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Abstract 

In this note, the statistical mechanics of SOS (solid-on-solid) 1-dimen- 
sional models under the global constraint of having a specified area between 
the interface and the horizontal axis, is studied. We prove the existence of 
the thermodynamic limits and the equivalence of the corresponding statistical 
mechanics. This gives a simple alternative microscopic proof of the validity 
of the Wulff construction for such models, first established in Ref. 5. 
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We consider the SOS model defined as follows: To each site i of the 
lattice an integer variable hi is assigned which represents the height of 
the interface at this site. The energy i?jv({/i}) of a configuration {h} — 
{ho,hi, ...,hN}, in the box < i < N, oi length N, is equal to the length of 
the corresponding interface 

N 

HNm) = 12i^ + \hi-hi.i\) (1) 

i=l 

Its weight, at the inverse temperature /3, is proportional to the Boltzmann 
factor exp[-pH{{h})]. 

We introduce the Gibbs ensemble which consits of all configurations, in 
the box of length N, with specified boundary conditions ho — and h^ — Y- 
The associated partition function is given by 

Z^{N,Y) = ^e-«'^»5(M5(/iiv (2) 
{h} 

where the sum runs over all configurations in the box and 6{t) is the dicrete 
Dirac delta = 1 if t = and 6{t) = otherwise). We define the 

corresponding free energy per site as the limit 

where y — —tanO, the slope of the interface, is a real number. This free 
energy is called the projected surface tension. The surface tension, which 
represents the interfacial free energy per unit length of the mean interface, is 

t(^) =cos^ Tp(-tan^) (4) 

We introduce also a second Gibbs ensemble, which is conjugate to the 
previous ensemble, and whose partition function, in the box of length A^, is 
given by 

Z^{N, x) = Y^ e-^^«'^»e'^^'^^5(/io) (5) 
{h} 
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where x & IR replaces as a thermodynamic parameter the slope y and ho — 0. 
We define the associated free energy as 

Theorem 1. Limits (3) and (6), which define the above free energies, 
exist. The first, Tp, is a positive bounded even convex function of y. The 
second, (p, is a bounded above even concave function of x. Moreover, Tp and 
— </? are conjugate convex functions, i.e., they are related by the Legendre 
transformations 

- (fix) = sup [xy - Tp{y)] 
y 

Tp{y) = sup [xy + cp{x)] (7) 

X 

Proof: The validity of the above statements is well known. See for in- 
stance Refs. 2, 3 for a proof of these results in a more general setting. 

The convexity of Tp is equivalent to the fact that the surface tension 
r satisfies a stability condition called the triangular inequality (see Refs. 
2, 3). Relations (7) between the free energies express the thermodynamic 
equivalence of the two ensembles (2) and (4). These relations imply that 
the curve z = (p{x) gives, according to the Wulff construction, or its modern 
equivalent the Andreev construction, the equilibrium shape of the crystal 
associated to our system. 

The function ip{x) defined by (6) is easily computed by summing a geo- 
metrical series. One introduces the difference variables 

Hi = hi_i - hi (8) 

ioT i — 1, N, so that the partition function factorizes and one obtains 

ip{x) = l-p-Hn ^ e-'^l"l+^^" (9) 
ne 'Ml 

The explicit form of this function is 

( ) 1 1 sinh^ 

^ ^ cosh P — cosh Px 
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if — 1 < X < 1 and (f{x) — — oo otherwise. 

We next define two new Gibbs ensembles for the system under consid- 
eration. In the first of these ensembles we consider the configurations such 
that = 0, which have a specified height at the origin Hq = M and which 
have a specified volume V between the interface and the horizontal axis, this 
volume being counted negatively for negative heights: 

N 

v^v{{h})=j:h, (11) 

i=0 

The corresponding partition function is 

Zs{N, V,M)=Y: e-'3^(W)<5(/i;v)5(V^({/i}) - V)6{ho - M) (12) 
{h} 

The second of these ensembles is the conjugate ensemble of the first one. Its 
partition function is given by 

Z,{N, u,fi) = J2 e-^^«''»e'^"(^«'^»/^)+'^'^'^°(5(/i^) (13) 
{h} 

where u & M and n E M are the conjugate variables. Our next step will be 
to prove the existence of the thermodynamic limit for these ensembles and 
their equivalence in this limit. 

Theorem 2. The following limits exist 

ip3(v,m) ^ lim --^ In Z3(N,vN^,mN) (14) 

ip4{u,ii) = lim --^ \nZ4{N,u,iJ,) (15) 

They define the free energies per site associated to the considered ensem- 
bles as, respectively, convex and concave functions of the v and u variables. 
Moreover, ip^ and are conjugate convex functions: 

— ip4{u, ji) = sup [uv + /im — ip3{v, m)] 

v,m 

ip3{v, m) = sup [uv + iim + ip4{u, /x)] (16) 



Proof: The crucial observation is the subadditivity property given in 
Lemma 1 below. Then we addapt known arguments in the theory of the 
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thermodynamic limit (see Refs. 4, 5). A more detailed discussion is given in 
the Appendix. 

Lemma 1. The partition function satisfy the suhadditivity property 
Z^{2N, 2{V' + V"),M' + M") 

> Zs{N, V\ M') Zs{N, V", M") e-^/^l^'^/P^-i) (17) 

Proof: In order to prove this property we associate a configuration {h} 
of the first system in the box of length 2iV, to a pair of configurations {/i} 
and {h}" of the system in a box of length N, as follows 

h2i ^h/i + h1 , i = 0,...,N 
h2i-i^K_^ + h'l , 1 = 1,. ..,N (18) 

Then = + ^jv = 0, /lo = /^o + K = M' + M" and 

N N N 

= 2\v{{h}')~+V{{hY)]-M''. 

This shows that the configuration {h} belongs to Z-i{2N ,2{y' + V") + 
M", M' + M"). Since H2N{{h}) = HN{{hy) + HN{{h}"), because = n'^ 
and n2i-i = n'l for i = 1, A'" — 1, as follows from (18), we get 

Z3{N, V, M') Z^{N, V", M") < Z^{2N, 2{V' + V") + M" , M' + M"). 

Then we use the change of variables hi = hi — [M" /{2N — 1)] for i = 
1, . . . , 2A" — 1, /lo = ^0; h2N — h2N — which gives 

Zs{2N, V + M", M) < e2/3|^"l/(2iv-i) ^^^gA^, V, M) 

to conclude the proof. 

Theorem 3. The functions ■03 and ■04 can he expressed in terms of the 
functions </? and Tp as follows 

1 

4>4{u, fj) = — (p{x + fj,)dx (19) 

u Jo 

Mv,m) = — Tp{ip'{x))dx (20) 

1*0 "'MO 
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where uq and satisfy 

1 /""o 1 

ip{x + iio)dx ifino + Uo) ^ V (21) 

Uq Jo uq 

— b()"o) - Vil^o + Uo)] = m (22) 

Uo 

Proof: We consider again the difference variables (8) and observe that 

N N 

V{{h})^J2hi^j:irii 

i=0 1=1 

and therefore 

N 

Z^{N,u,fj,)^]l[ e-^l"'l+^("/^)^"'+'^'^"') 
^=1 me'E 

Taking expression (9) into account it follows 

Z4{N, -u, ^) = exp - /3 ^ ip{—i + /i)) 

i=l -'^ 

and 

1 ^ u 1 ^ u u 

Mu, /^) = Jim - ^{-^ + ^^)^ Jim - ^ ^¥'(^i + /^) 
1=1 j=i 

which implies (19) in the Theorem. 

The function ips is determined by the Legendre transform (16). The 
supremum over u,^ is obtained for the value Uo, Ho for which the partial 
derivatives of the right hand side arc zero: v + {d/du)i{j4{uo, /J^o) = 0, v + 
{d/dn)ilj4:{uo, Ho) = 0. That is, for Uo, Ho which satisfy (21). 

Then, from (16), (19) and (21), we get 

ip3{v, m) = 2il)4{uo, Ho) - —[{Ho + Uo)(f{Ho + Uo) - HofiHo)] (23) 

1*0 
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Figure 1. Graphical interpretation of Theorem 3. 



The right hand side of (23) represents twice the area of the sector OBC in 
Fig. 1 divided by uq. But, it is a known property in the Wulff construction, 
that twice this area is equal to the integral in (20). Indeed, by using the 
relation (7) in the form 

(fi{x) = X(fi'{x) + Tp{(fi'{x)) 

in (20) and integrating by parts x(p'{x), we get 
2ipi{uo, no) = — Tp{ip {x))dx H [(/io + uo)ip{(Xo + uq) - Hovino)] 

Uq JiXQ Uq 

which together with (23) imphes the expression (19) in the Theorem. 

To interpret these relations, let us observe that the right hand side of (21) 

2 

represents the area ABC, in Fig. 1, divided by AC . Therefore, the values 
Uq and Uq, which solve (21) and (22), arc obtained hen this area is equal to 
V, with the condition, coming from (22), that the slope AB/AC is equal to 
m. Then, according to (19), the free energy %lj^{y,m) is equal to the integral 
of the surface tension along the arc BC, of the curve z = (j){x), divided by 
the same scahng factor AC — Uq. 

We conclude that, for large N , the configurations of the SOS model, with 
a prescribed area vN"^, follow a well defined mean profile, the macroscopic 



8 



S. Miracle-Sole and J. Ruiz 



profile given by the Wulff construction, with very small fluctuations. This 
follows from the fact that the probability of the configurations which deviate 
macroscopically from the mean profile is zero in the thermodynamic limit. 
The free energy associated to the configurations which satisfy the above 
conditions, and moreover, are constrained to pass through a given point not 
belonging to the mean profile, can be computed with the help of Theorem 
3. The corresponding probabilities decay exponentially as N ^ oo, as a 
consequence of the usual large deviations theory in statistical mechanics (see 
Ref. 6). 
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Appendix 

We give here a more detailed discussion of the proof of Theorem 2. We 
define: 

fn{v,m) = -^\nZs{2^,2'^v,2^m). 

For V and m of the form 2^^p, the subadditivity property with N — 2", 
V — V" — vN'^ and M' — M" — mN imphes that is a decreasing 
sequence: /„+i(v,m) < fn{v,m). Since this sequence is bounded from below 
its limit exits when n tends to infinity. Indeed 

{h} 

The R.H.S. of the above expression is easily computed by inroducing the dif- 
fernce variables (8) and we get that is bounded from below by ln(l — 

e~^). Let us notice that one can obtain a lower bound to Z3{N,vN'^,mN) 
by restricting the summation over the configuration such that (A^ — l)hi = 
vN^ — {l/2)mN for i — 1, . . . , N — 1. This gives that /„ is bounded from 
above by |v — m| + It;!. 

The existence of the limit for general N follows from standard argument 
in the theory of the thermodynamic limit (cf . 4) and we have: 



ipsiv, m) = inf 



^ lnZ2{N,vN\mN) 



(Al) 



To prove that ^3 is convex, we notice that the subadditivity inequality 
(17) with N = 2", V = viN^, V" = v-iN'^ and M' = M" = mN"^ gives: 

i^ii^ivi + V2,m)) < ]^i)i{vi,m) + ]^tl)^{v2,m) 
which applied iteratively implies: 

V^a^avi + (1 - a)v2-,m}j < aip3{vi,m) + (1 - a)ilj3{v2,m) 
for a of the form 2~*p and < a < 1. For such a we obtain analogously 

ilj3(v,ami + (1 - a)m2^ < ail)3{v,mi) + (1 - a)i})3{v,m2) 

by applying the subadditivity inequality (17) with N — 2", M' — rriiN^, 
M" = m2N^ and V = V" = vN'^ 
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Since V's is bounded, it follows that for all m, ips can be extended to a 
convex Lipshitz continuous function of the real variable v. 

To prove the the existence of the limit (15) and relations (16), we intro- 
duce 

ZtiN,u,fi)= sup [eWW^M ^^(^^y^^y 
v,Me S 

and proceed, as in the Appendix of Ref. 5, to study the thermodynamic limit 
for this quantity. Let 

■04 (w, /i) — sup [uv + iim — ip3{v, m)]. 

v,m 

According to (A.l), we have 
for all y e so that 

Z+{N,u,n) < e^^^^H.a (A2) 

On the other hand for any 6 > and sufficiently large N one can find, 
V = vN"^ and M = mN such that 

> ^pN[n(u,n)-5\ 

and therefore 

Z+{N,u,fj) > eW:KM)-<5]. (^.3) 
Inequalities (A. 2) and (A. 3) imply that 

^^Z+{N,u,pl) = -r^{u,pL). (A4) 

N-^oo ply 

We shall now prove that the thermodynamic limit (15) exits and gives the 
same quantity following the argument of Theorem 2 in Ref. 5. First, we notice 
that 

y,Me 

which implies 

Z+{N,u,pi)<Z^{N,u,pL) (A5) 
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Moreover, the inequality 

used with u — u' and u — u", p, — /i' and p, — /i" gives for any u" < u < u' 
and any /i" < /i < /i' : 

+ Z* Ey<o e/'(""-")(^/^) Em">o e'^^'^"^')'" 

+ Z* Ey>o e/^(«""')(^/^) Em<"o e^^i^"-^)^ {A.6) 

+ Z* Ey<o e/^(""-")(^/^) Em"<o e/^^'^""'^)^ 



(l_e-'9(^''/'^))(l-e-'3(AM)) 

where 

= sup [Z4+(iV, m', /i'), ^4^ Ai"), ^4+(^, /^'), ^4^(^, U", /^")] 

< Ali = min[ii' — u,u — u"] 
< Au — min[//' — /i, /i — /i"] 

By referring to (A. 4), and to the continuity of ■04, the inequahties (A. 5) and 
(A.6) imply (15) and (16). 



